Multiplexing of discrete chaotic signals in presence of noise 
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In this paper, multiplexing of discrete chaotic signals in the presence of noise is investigated. 
Existing methods are based on chaotic synchronization which is susceptible to noise and parameter 
mismatch. Furthermore, these methods fail for multiplexing more than two discrete chaotic signals. 
We propose two novel methods to multiplex multiple discrete chaotic signals based on the principle of 
symbolic sequence invariance in the presence of noise and finite precision implementation of finding 
the initial condition of an arbitrarily long symbolic sequence of a chaotic map. 
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I. MOTIVATION 



In this paper, we investigate multiplexing of chaotic 
signals in the presence of noise. The fact that chaotic 
signals have a lot of redundancy in them is exploited 
here. We first motivate the application and then provide 
new methods. 

Multiplexing of signals is a very important requirement 
in most communication systems. Consider the scenario 
where there are multiple signals from multiple senders 
to be transmitted to multiple receivers, but there exists 
only one communication channel that can transmit only 
one signal at any given time. In such a scenario, it would 
be beneficial if all the signals are "added" in a special 
way to create a single composite signal for transmission 
across the communication channel. This single composite 
signal is "separated" to the respective signals in a lossless 
fashion at the other end of the channel. Note that there 
is also noise which is invariably added at the channel. 
This scenario can also occur in transmission of neuronal 
signals from different parts of the brain to various parts 
of the body through a single channel. 

For linear communication systems, standard ways such 
as frequency division multiplexing (different signals are 
allocated different parts of the frequency spectrum) and 
time division multiplexing (different signals are allocated 
different time slots for transmission) are used to increase 
the information capacity of the channel [1]. Nonlinear 
chaotic oscillators are increasingly being used in com- 
munications since it offers a potential advantage over 
conventional classical methods in terms of noise perfor- 
mance [2]. Multi-user chaotic communications has be- 
come a hot topic of research in recent times [3] . It is also 
potentially useful in spectrum-spreading communication 
systems. Hence there is a need for multiplexing chaotic 
signals. 
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II. EXISTING WORK AND THEIR 
LIMITATIONS 

There has been some work in multiplexing chaotic sig- 
nals. For the first time in 1996, multiplexing of chaos 
using chaotic synchronization was investigated in a sim- 
ple map and an electronic circuit model by Tsimring and 
Sushchik [4]. Liu and Davis [5] used a scalar signal to 
simultaneously synchronize two different pairs of chaotic 
oscillators. They called this method as dual synchro- 
nization. However, there are several limitations of this 
method. They derive a condition for dual synchroniza- 
tion which holds only for certain discrete chaotic signals 
(maps) and for certain values of the coupling coefficients. 
The notable omission is the binary map (Bernoulli shift). 
They show that the binary map does not satisfy the 
condition for dual synchronization for any value of the 
coupling coefficients. Thus, chaotic signals from the bi- 
nary map can't be multiplexed by their method. Another 
limitation is that their method can only work with two 
chaotic signals. It is not known whether the method can 
be extended to multiple signals (more than two) from 
different maps. 

There has been more work in multiplexing chaotic sig- 
nals from continuous chaotic systems (flows). Liu and 
Davis [5] extend their work to multiplex signals from 
delay-differential equations. Further progress has been 
made by Di Ning et. al. [6] who extend Liu's method for 
3D continuous chaotic systems (Lorenz and Rossler sys- 
tems). This has been further improved by Salarich and 
Shahrokhi [7] who make use of a time-varying output 
feedback strategy to achieve dual synchronization with- 
out the need for all the master states (only a linear com- 
bination of the master states are enough). 

All the developments were only for two chaotic systems 
until the work of Salarieh and Shahrokhi [8] who suc- 
ceeded in multiplexing more than two continuous chaotic 
signals using chaotic synchronization via output feedback 
strategy. They derive a necessary condition for multi- 
synchronization and demonstrate the algorithm to the 
Chen-Lorenz- Rossler and the Duffing- Van der Pol contin- 
uous time chaotic dynamical systems. However, none of 
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these methods work for multiple discrete chaotic signals. 
Another serious limitation is that chaotic synchroniza- 
tion is susceptible to noise and to parameter mismatch. 
Even one percent of parameter mismatch leads to 8% of 
synchronization error and one percent of noise results in 
a synchronization error of 4% as reported by [5]. 

Vaidya's method [9] which multiplexes more than two 
discrete chaotic signals remains as the latest development 
on multiplexing chaotic signals from discrete chaotic sig- 
nals (ID maps). Vaidya's method does not use chaotic 
synchronization and is fundamentally different from the 
previous approaches. We will describe the method and 
its drawbacks in Section IIIC since we are going to use 
some of the ideas from this method to improve upon it. 
In principle, it is possible to extend the methods that is 
proposed in this paper to flows and to higher dimensional 
chaotic dynamical systems, but these will not be pursued 
here. Since Liu and Davis' method does not work for the 
standard binary map, we shall consider chaotic signals 
from the standard binary map with randomly chosen ini- 
tial conditions. 



A. New Approach 

Our approach considers two different scenarios as 
shown in Figure 1. In Scenario 1, the communicating 
channel is noisy and there is no control on the noise that 
is added during transmission of the signal (Figure 1(a)). 
However, it shall be assumed that the magnitude of the 
noise is limited (we shall give conditions on the magni- 
tude of noise that is allowed by our methods) . The noise 
is uniformly distributed. The signals are chaotic and the 
noise that is added during transmission is uncontrolled 
but limited in magnitude. 

In Scenario 2 (Figure 1(b)), the communication chan- 
nel is lossless, but the noise is added at the sender. The 
noise is assumed to be of the same magnitude as the 
chaotic signals and uniformly distributed. But the way 
the noise is added is under control. This scenario cor- 
responds to steganography [17] or cryptographic appli- 
cations where the noise could be the "payload" to be 
secretly transmitted. 

Recently, Vaidya [9] suggested a novel multiplexing al- 
gorithm for ID discrete chaotic signals in the presence of 
noise. We shall call this as Method 1 and review it briefly 
and list some of its limitations. A new method (Method 
2) will be proposed which overcomes some of the limita- 
tions of Method 1. Both these methods are solutions to 
Scenario 1. A novel method (Method 3) is proposed for 
Scenario 2. 



III. METHOD 1: VAIDYA'S NOISE RESISTANT 
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(a) Scenario 1: Channel is lossy, noise is additive and limited 
in magnitude (Methods 1 and 2). 
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(b) Scenario 2: Channel is lossless, but noise which has the 
same magnitude as the signal, is added in a special way at 
the sender (Method 3). 

FIG. 1: Multiplexing of chaotic signals in the presence of 
noise. 



does not make use of chaotic synchronization like the pre- 
vious approaches. Vaidya proposes a noise-resistant ver- 
sion of the Tent map. Since we are going to deal mainly 
with the standard binary map, a minor modification leads 
to a noise-resistant version of the binary map. It is given 
by the following set of equations: 



2x, 

2x-hg, 
0, 



0<.<| 



< X < p 



p < X < p + q 



ip 



2x~2{p + q), p + q<x<— + q 
2x-{2p + q), y -t-g < a; < 2p + 9 



= 0, 



2p + q < X < I. 



Figure 2 depicts the noise resistant binary map (de- 
noted by Tnoiseres)- Vaidya establishes that there ex- 
ists a conjugacy between the ordinary binary map and 
Tnoiseres- Givcu any chaotic signal on the binary map 
(a chaotic signal is trajectory on the map for a given 
initial condition), one can find the equivalent signal on 
Tnoiseres- P (0 < p < 0.5) and <z (0 < (? < 0.5) can be 
chosen such that 2p+2q= 1. 

The symbolic sequence is defined as follows: 



Vaidya's method [9] is a solution to multiplexing of 
chaotic signals in the presence of noise (Scenario 1) which 



S{xi) = 0, < Xi < p q 
= 1, p + q < Xi < I. 
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FIG. 2: Tnoiseres- Noise-resistaiit Binary Map. 

Here, X = {xiY~J!^ is the chaotic trajectory (or chaotic 
signal) starting from an initial condition xi. S{X) de- 
notes the symbolic sequence for the entire trajectory. 

A. Noise-resistance 

For any given chaotic signal X on Tnoiseres, if noise 
N = {niY~^™ is added such that each rii satisfies < 
rii < q, then it can be seen that the symbolic sequence 
remains unchanged: 

S{X) = Six + N). (1) 

The signal X is transmitted at the sender and the re- 
sulting signal Z = X + N is what is received at the 
receiver. However, because of the above property of sym- 
bolic sequence invariance, we can compute S{X) {=S{Z)) 
and iterating backwards on the map, we can find the ini- 
tial condition xi. Knowing xi, we can easily compute 
X2, ■ ■ ■ , Xm and thus we can recover X losslessly. Fur- 
thermore, we can also compute N = Z X accurately. 
Thus we have recovered both the original chaotic signal 
and the noise. This noise resistance property is provided 
by a non-zero value of q. The larger the q, the higher is 
the resistance to noise but at the same time the length of 
the signal has to be longer in order to determine the ini- 
tial condition xi accurately from the symbolic sequence. 

B. Cascading Noise-resistant Maps 



sequence invariance is maintained at each step of cascad- 
ing. Thus, the symbolic sequence of Z is used to decode 
Xi and Ni where Ni is the sum of {X2, X3, . . . ,Xk} and 
N. The symbolic sequence of Ni is the same as that 
of X2 and hence X2 can be decoded. This procedure 
is repeated until all the signals are losslessly recovered 
along with N. Vaidya successfully applies this method 
to multiplex 100 chaotic signals. 

C. Drawbacks of Method 1 

The drawbacks of Method 1 are as follows: 

1. Given chaotic signals from various chaotic maps, 
one has to find the corresponding signals in the 
noise- resistant binary/tent map using topological 
conjugacy. This is cumbersome and although we 
can have a finite precision implementation of the 
chaotic maps and the noise-resistant binary map, 
it is difficult to develop finite precision implemen- 
tation of the topological conjugacy. For example, if 
the signal X is from the logistic map, then the topo- 
logical conjugacy will involve trigonometric func- 
tions which would have to be expressed as trun- 
cated infinite series. This may lead to errors in the 
recovery of X. 

2. The amount of noise that can be added reduces 
exponentially as the number of chaotic signals to 
be multiplexed increases linearly. 

3. It is in principle possible to extend the idea of noise 
resistant maps to other chaotic maps like the logis- 
tic map. However, for each new map, the equations 
have to be worked out explicitly. 

We are motivated to invent new methods of multiplex- 
ing which will circumvent the above problems. While ex- 
ploiting the idea of symbolic sequence invariance under 
the addition of noise, we would like to devise a method 
which will work for any ID chaotic unimodal map (and 
generalizable to other kinds of maps and higher dimen- 
sional ones) without the necessity of topological conju- 
gacy. The scenario where the magnitude of noise is equal 
to that of the signal also needs to be addressed. 



Vaidya goes one step further and defines a cascade of 
such noise-resistant maps. To add another chaotic sig- 
nal Y to X, Vaidya defines a similar noise-resistant map 
which maps [0, q) onto itself. It is self-similar to Tnoiseres- 
Thus, he defines a whole cascade of noise-resistant maps, 
all of which are self similar to the original one. The 
domain of succeeding maps reduces exponentially. For 
further details, please refer to [9] . 

With these cascade of maps, one could now add a whole 
family of chaotic signals {Xi,X2, . . . ^Xk} to one noise 
signal TV (on the channel) with magnitude dictated by 
the number of maps, to yield the signal Z . The symbolic 



IV. METHOD 2 

The key idea of Method 1 is the notion of symbolic 
sequence invariance. As long as we ensure that the sym- 
bolic sequence of the original chaotic signal X is unaf- 
fected by adding the noise N (uniformly distributed), 
the resulting signal Z ^ X + N has the same symbolic 
sequence as X {S{Z) = S{X)). Then, given this arbitrar- 
ily long symbohc sequence S{X), the problem reduces to 
determining its initial condition and iterating this initial 
condition to obtain the entire chaotic signal X . Once X 
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is determined, one could subtract X from Z to obtain 
the noise signal N. 

In order to find the initial condition from an arbitrarily 
long symbolic sequence of a chaotic map, we make use of 
GLS-coding. Generalized Luroth Series or GLS-coding 
for short, is a new entropy coding algorithm [10] that 
achieves the Shannon's entropy rate for lossless data com- 
pression. The idea of GLS-coding is to first embed the 
stochastic (binary) i.i.d source into the appropriate GLS 
and then treating the message as a symbolic sequence, 
the initial condition is determined by a backward itera- 
tion. A finite precision implementation of GLS-coding is 
described in the Appendix. Using such an implementa- 
tion, it is possible to determine accurately the initial con- 
dition of an arbitrarily long symbolic sequence. Our im- 
plementation is for the skew-tent and skew-binary maps 
and can be extended to other maps. The standard tent 
map and binary map are part of this family. 

The algorithm for Method 2 is described as follows: 



9. Once we have Dnoisy = D, we can recover the 
symbolic sequences of each of the k chaotic sig- 
nals and thereby recover the k initial conditions 
{Xi(l),X2(l), . . .,Xk{l)} by GLS-coding (see Ap- 
pendix). 



10. From the initial conditions, the k chaotic signals 
can be recovered. 



1 1 . The noise signal can be recovered at the receiver by 
computing Z from D and computing N — Znoisy — 
Z. 



0' 
0.25 



k=1 



0.75 



1. Let Xi,X2, ■ ■ ■ Xk be k chaotic signals of length m 
to be multiplexed. Each of these signals are ob- 
tained from distinct initial conditions (randomly 
chosen) on the standard binary map. 



2. Compute the symbolic sequence 
{S{X,il)),S{X,i2)),...S{X,{m))}l 
tion S(.) is defined as follows: 



{sx^l 



The func- 



S{x) =0, < a; < 0.5 
= 1, 0.5 < a; < 1. 



(2) 
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FIG. 3: Method 2: For the cases = 1 and k = 2. The points 
marked 'X' are transmitted depending on the symbolic se- 
quences of the chaotic signals. Owing to noise at the channel, 
the received signal will also be uniformly distributed in the 
range [0, 1). This explains Figure 5(c). 



3. Compute {S'i}*=f ^ {< SXi{l)SX2{l). . . 
SXk{l)>2, < SXi{2)SX2{2)...SXk{2) >2 
... < SXi{m)SX2im)...SXk{m) >2}. Here, 
< . >B denotes a number that is base-B represen- 
tation. 

4. Compute D = {A}^^?" where D, S^ >io. 

5. Compute Z = {Z,}™ where Z, = ^§^. 

6. Transmit Z across the channel. 

7. Receive Z^ozsy ^ Z + N. Here = {N,}™ where 
each Ni is uniformly distributed noise in the range 
^_2-(fe-Hi)^_,_2-(fc+i))^ 

8. At receiver, compute Dnoisy = Y^^Znoisy\ where 
[(.)J is the floor operation which computes the 
maximum integer that is less than the argument. 



Note that Dnoisy = D. This is because D < 
2''Znoisy < {D + 1). The fioor operator makes this 
equal to D since D is always a positive integer. This 
is where we have made use of the fact that the sym- 
bolic sequence is invariant in spite of noise. Here 
D has the information of the symbolic sequence of 
all the k chaotic signals. 
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FIG. 4: Method 2: Multiplexing of fc = 10 chaotic signals 
in the presence of one noise signal TV - (a) Top left: The 7th 
Chaotic Signal X7{.). m — 1000, but only first 100 shown, (b) 
Top right: Phase Portrait of X7{.), (c) Bottom: Histogram of 
Xr(.). 
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FIG. 5: Method 2: Multiplexing of chaotic signals in the 
presence of noise -(a) Top left: Z signal (showing first 100 
values only), (b) Top right: Phase portrait of noise signal AT, 
(c) Bottom left: Phase portrait of Znoisy, (d) Bottom right: 
Histogram of Z„oisy 



magnitude of noise (2~'^) that can be tolerated at the 
channel, goes down exponentially. Method 3 overcomes 
this limitation. The noise magnitude can be equal to the 
signal. However, we can no longer operate in Scenario 1. 
We assume that we have control on the "way" the noise is 
added (noise is still assumed to be uniformly distributed) 
and that the channel is lossless (Scenario 2). 
Method 3 is described as follows: 

1. Let Xi,X2, ■ ■ ■ Xk be k chaotic signals of length m 
to be multiplexed. Each of these signals are ob- 
tained from distinct initial conditions on the stan- 
dard binary map. 

2. Let noise signal be = where each Ni is 
independent and identically distributed (uniform) 
in the range (0, 1). Noise signal N is independently 
generated but available at the receiver. 

3. Given two signals A = {A(i)}*=™ and B = 
{B{i)}l'z^ where ^ is a chaotic signal {B can be 
anything), we define the operation A + B.S{A) = 
{A{i) + B{i).S{A{i))}lZT as follows: 



A. Experimental Simulations 

As a simple example. Figure 3 shows the points trans- 
mitted for the case k = 1 and k = 2. 

Method 2 was experimentally simulated for /c = 10 
chaotic signals of length m = 1000 each. They were 
all generated by randomly chosen initial conditions on 
the standard binary map. Figure 4 shows the seventh 
chaotic signal Xj. The phase portrait and the histogram 
are also shown. Figure 5 shows Z, the phase portrait of 
noise signal A'^, phase portrait of Znoisy after the addition 
of noise and the histogram of Znoisy AH the 10 chaotic 
signals and the noise signal were successfully recovered in 
a lossless fashion at the receiver. We have used the finite 
precision implementation of GLS-coding as proposed in 
Appendix. This confirms the efHcacy of Method 2. 



V. METHOD 3 

The biggest advantage of Method 2 is that in princi- 
ple it works for any dynamical system. As long as we 
know the Markov partitions of the dynamical system, we 
can define the symbolic sequence and hence use Method 
2. There is no need of using topological conjugacy or 
construction of special noise-resistant maps like Method 
1. However, one needs to develop an analog of GLS- 
coding (i.e. finding initial condition for an arbitrary long 
symbolic sequence of the dynamical system as given in 
Algorithm 1 in Appendix) for the method to work. 

Method 2 works in the presence of a lossy channel, the 
noise being additive, but the magnitude of noise that is 
tolerated depends on the number of signals being mul- 
tiplexed. As the number of signals (fc) increases, the 



A{i)+B{i).S{A{i)) = A{i)-B{i), iiS{A{i))=0 
= Aii) + Bii), i{S{A{i)) = l. 

where S{.) is defined in Equation 2. 

4. Compute the following signals: 

Xi+N.S{Xi) + l 



Zk 

z 



X2 + Zt.S{X2) + 1 

3 

X3 + ^2.5(^3) + 1 



Xk + Zk-i.S{Xk) + 1 
3 

Zk- 



5. Transmit Z on the lossless channel. Note that the 
dynamic range of Z is [0, 1). 

6. Receiver receives Z. By symbolic sequence invari- 
ance, we have the following identities: 



S{Z) = S{Zk) 
S{Zk-i) 



S{Xk) 
S{Xk-i) 



S{Z,) = S{X,). 

7. We start with Z and compute S{Zk). By the first 

identity, we have S{Xk)- GLS-coding is applied to 
determine ^fe(l). Hence Xk is recovered losslessly. 
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Knowing Xk and S{Xk)^ we can compute Zk-\ by 
the following equations: 



Zk_^(i) = ^ 3Z(«) + 1, if 5(Xfe(z.)) = 

= 3Z(i) - Xfe(0 - 1, if 5(Xfe(0) = 1. 

8. Knowing 2fe-i) we repeat the procedure to extract 
and Zk-1- This is repeated until we have 
extracted all the chaotic signals and the noise signal 
A'^. Note that noise can be thought of as Zq and 
the same procedure applies. 




FIG. 6: Method 3: Multiplexing of k = 24 chaotic signals 
in the presence of one noise signal - (a) Top left: The 17th 
Chaotic Signal Xi7{.). The length m — 1000, but only the 
first 100 shown, (b) Top right: Phase Portrait of Xi7{.), (c) 
Bottom left: Noise signal A''(.), (d) Bottom right: Phase 
Portrait of N{.). Note that magnitude of N is the same as 
that of X. 



A. Experimental Simulations 

Method 3 was experimentally simulated for fc = 24 

chaotic signals of length m = 1000 each and one noise 
signal of the same length. The chaotic signals were all 
generated by randomly chosen initial conditions on the 
standard binary map. Figure 6(a) shows the 17th chaotic 
signal Xir. The phase portrait is shown in Figure 6(b). 
Figure 6(c) shows the noise signal N which has the same 
magnitude as that of the chaotic signal. The phase por- 
trait of noise signal N is shown in Figure 6(d). The final 
signal that is transmitted on the lossless channel Z is 
shown in Figure 7(a). Its phase portrait and histogram 
are shown in Figure 7(b) and (c) respectively. Note that 
in this method, the final signal Z docs not have uniform 
distribution although the individual chaotic signals and 



the noise signals arc uniform. This is because of the spe- 
cial way in which the signals were added. 

All the 24 chaotic signals and the noise signal were 
successfully recovered in a lossless fashion at the receiver. 
We have used the finite precision implementation of GLS- 
coding as proposed in Appendix. This confirms the effi- 
cacy of Method 3. 




z 



FIG. 7: Method 3: Multiplexing of chaotic signals in the pres- 
ence of noise - (a) Top left: Z signal generated by Method 3 
is not uniformly distributed. This is transmitted on the loss- 
less channel, (b) Top right: Phase portrait of Z, (c) Bottom: 
Histogram of Z showing that it is bimodal. 



VI. REMARKS ON THE THREE METHODS 

The following observations can be made on the three 
methods: 

1. The idea of symbolic sequence invariance is the key 

to the success of all the three methods. The way 
this idea is implemented is different in the three 
methods. 

2. The way noise is handled is the same in Methods 

1 and 2 since there is no control on the noise in 
Scenario 1. Scenario 2 is much more restrictive in 
terms of noise. 

3. All the three methods rely on the finite precision 
implementation of GLS-coding, i.e. finding the ini- 
tial condition given an arbitrarily long symbolic se- 
quence of the chaotic signal (refer to Appendix). 
The method of extracting the symbolic sequence 
from a "noisy non-linear system" and performing 
GLS-coding is analogous to filtering out noise in 
linear systems by means of integration or other lin- 
ear filters (for eg., low pass filters). 
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4. Methods 2 and 3 can be easily extended to tent 
map, skew-tent map, logistic map and other uni- 
modal maps. It is also possible to extend the meth- 
ods to non-unimodal ID maps and possibly higher 
dimensional maps. The key is to find an analog 
of GLS-coding in those cases, which we believe is 
possible. 

5. Method 1 has a hint of using the idea of 'forbid- 
den symbol' [11-14] by allocating the length q on 
the interval which is never used by the map. This 
method can be potentially used for error correction 
and detection. 

6. In Method 2, just by observing the signal on the 
channel, no information can be gleaned. The distri- 
bution is uniform and the phase portrait is also ran- 
dom looking. The fact that multiple chaotic signals 
have been embedded is not obvious. This property 
enables it to be used in steganography or informa- 
tion hiding. In LSB-steganography, the Least Sig- 
nificant Bit (LSB) of natural signals is replaced by 
the secret (noise or noise- like). Method 2 is doing 
the reverse: MSB-steganography, where the Most 
Significant Bit (MSB) of the secret (noise or noise- 
like) is replaced by the symbolic sequence of the 
chaotic signal. 

7. Method 3 can handle any number of chaotic signals 
and one noise signal. However, in practice there 
will be limitation on the number of signals owing 
to finite precision since we are rescaling the range 
of the signals to [0, 1) (by addition of 1 and division 
by 3). 

8. The methods show that chaotic signals are highly 
redundant and hence robust to noise. As long as 
the symbolic sequence is preserved, the actual sig- 
nal can be distorted to a great deal. Also, forward 
iteration of chaotic dynamical systems exhibits sen- 
sitive dependence on initial conditions, but back- 
ward iteration shows resistance to noise. These fea- 
tures are not exhibited by random/stochastic sig- 
nals. This probably makes a strong case for why 
biological systems may use chaotic signals for trans- 
mission of information. Neuronal signals may use 
similar mechanism for robust transport of informa- 
tion. 

9. The above methods will not work for purely ran- 
dom signals or for non-chaotic signals since there 
is no way we can construct the entire trajectory by 
knowing the symbolic sequence. The redundancy 
of chaotic signals is necessary. At the same time, 
chaotic signals appear "random" in distribution. 

10. In contrast with earlier work on multiplexing dis- 
crete chaotic signals using chaotic synchronization. 



our methods can multiplex more than two discrete 
chaotic signals. There are no coupling coefficients 
used in our methods and hence there is no condi- 
tion to be satisfied for multiplexing. As noted pre- 
viously, earlier methods based on chaotic synchro- 
nization are vulnerable to noise and parameter mis- 
match. Our methods are completely robust to any 
amount of parameter mismatch since our methods 
do not rely on chaotic synchronization. Methods 1 
and 2 can tolerate noise but limited by the number 
of signals that is added. However, in Method 3, the 
noise magnitude is the same as that of the signals. 

11. The methods that we have developed can be poten- 
tially used in communication protocols, cryptogra- 
phy and steganography applications. 

12. There is no violation of Shannon's theorems for in- 
formation transmission in any of the methods. By 
transmitting the entire trajectory, we are sending 
lots of bits, much more than actually required for 
sending only the initial condition. These methods 
are not meant for compression of data. These are 
mechanisms to exploit the inherent redundancy in 
chaotic signals in spite of noise. 



VII. CONCLUSIONS AND OPEN PROBLEMS 

In this work, GLS-coding was used for multiplexing of 
chaotic signals in the presence of noise. By using the idea 
of symbolic sequence invariance, we were able to "add" 

several chaotic signals and "separate" them losslcssly at 
the receiver. We can either have a lossy channel but with 
limited noise (Methods 1 and 2) or have a lossless chan- 
nel with noise having the same magnitude as the signal 
"added" in a very special way at the sender (Method 3). 
An open problem is to investigate whether one can have 
both features in a single method. 

The inherent redundancy and structure in chaotic sig- 
nals which otherwise appear random in probability dis- 
tribution can be harnessed for robust communication of 
information. It is quite likely that such efficient mech- 
anisms (or similar ones) of handling noise in dynamical 
systems are already being employed in naturally occur- 
ring physical and biological systems. 

Compared to existing methods of multiplexing discrete 
chaotic signals, our methods arc significantly superior in 
all respects. The newly proposed methods can handle 
multiple signals from multiple maps (including Bernoulli 
shift or the binary map which was not possible by the 
method of Liu and Davis), completely robust to param- 
eter mismatch and good noise resistance capability. 
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Appendix: Finite Precision Implementation of 
GLS-coding 

The idea of GLS-coding [10] is to first embed the 
stochastic (binary) i.i.d source into the appropriate GLS 
and then treating the message as a symbolic sequence, 
the initial condition is determined by a backward itera- 
tion. This actually results in an interval (since the sym- 
bolic sequence is finite in length) and the mid-point of the 
interval is used as the initial condition. As the length of 
the symbolic sequence increases, the interval in which the 
initial condition is going to lie shrinks in size. This cre- 
ates problem in performing the backward iteration on a 
finite precision computer as the two ends of the interval 
come very close to each other and at some point it would 
be no longer possible to continue with the backward iter- 
ation. This problem needs to be addressed by some kind 
of re-normalization or re-scaling of the interval, in order 
for the method to be useful for encoding long sequences. 
Another problem with the method is that there is a long 
encoding delay. No output can be written/sent unless 
the entire initial condition is determined, which happens 
only after all the input bits are encoded. Luca's chaotic 
compression method [15] also has similar problems but is 
not addressed by them. 

There has been efforts to address both these problems 
for Arithmetic Coding [16]. Since GLS-coding is an ex- 
tension of Arithmetic Coding, these methods could be 
used. 

The idea is as follows. As soon as the interval com- 
pletely lies to the left of p (for the standard tent map 
and binary map, p ~ 0.5), the final initial condition will 
have a '0' in its binary expansion ('1' if the interval is 
completely to the right of p). Hence, this can be writ- 
ten as output and the current interval can be doubled 
in length. This ensures that the two ends of the inter- 
val will never come close to each other. At the end of 
all the iterations, the mid-point of the final interval is 
written as output. The only case in which this method 
would fail is when the interval straddles 0.5 at every it- 
eration. The probability of this happening exponentially 
decreases with each iteration. To handle this special case, 
there can be a check on the size of the interval and once 
it reduces to certain value, the encoder is forced to gen- 
erate an output and the interval is reset to [0,1). The 
iteration starts afresh for the next incoming bits. This 
would increase the size of the compressed file slightly as 
we are not encoding the entire symbolic sequence to de- 
termine a single initial condition, but the increase in size 
is negligible for long sequences. A similar technique was 
used in the IBM Q-coder [16] to handle this problem. 
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The algorithm for encoding is described in Algorithm 1. 
The decoding algorithm is similar and is omitted here. 

Algorithm 1 Finite precision GLS-coding 

1. Input: Binary message M of length A'^ from stochastic 
binary i.i.d source X. 

2. Compute probability of '0' as p = 

total number of zeros in M 

JV • 

3. Embed source X in to a GLS: construct GLS with par- 
titions [0,p) corresponding to symbol '0' and partition 
\p, 1) corresponding to symbol '1'. 

4. Initialize interval [L, U) to [0, 1). Initialize Tol = 10"*. 

5. Initialize k = \. 

6. Input the fc-th bit from M. 

7. If the bit is 0, then set: 

L ^ L/p 
U ^ U/p. 



For the multiplexing of chaotic discrete signals pro- 
posed in this paper, we do not need to compute p (in 
step 2 of Algorithm 1), but can directly set it to p = 0.5 
since we are using the standard binary map. The al- 
gorithm described here can be extended easily to find 
the initial condition from an arbitrarily long symbolic 
sequence of other chaotic maps (for eg., to find the ini- 
tial condition on the skew-tent map, equation in step 7 
needs to be modified appropriately) . A similar extension 
can be done for the Logistic map as well. 



else, set: 

L ^ (L-p)/{l-p) 
U ^ {U-p)/{l-p). 

Set fe ^ fe -I- 1. 

8. if < L, < 0.5, then Output bit '0' and set: 

L 2L 
U <- 2U. 

if 0.5 <L,U <\, then Output bit '1' and set: 

L <- 2L - 1 
U <- 2U-1. 

9. if (U — L) < Tol, Output Xic — ^"^^^-^ in binary repre- 
sentation and reset [L, U) — [0, 1). 

10. If fe < A'^, go to step 6, else continue to step 11. 

11. If L < 0.5 < U, then Xie = 0.5 else x^e = 

12. Output Xic in binary representation. 

13. Output p, N and mode used as overheard information. 
The precision of p can be chosen conveniently. 



